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Abstract 

We provide a method to correct the observed azimuthal anisotropy in heavy-ion collisions for the event plane res- 
olution in a wide centrality bin. This new procedure is especially useful for rare particles, such as Q baryons and 
J/t// mesons, which are difficult to measure in small intervals of centrality. Based on a Monte Carlo calculation with 
simulated V2 and multiplicity, we show that some of the commonly used methods have a bias of up to 15%. 
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1. Introduction 

Azimuthal anisotropy is one of the key observables 
to study the properties of matter created in high energy 
heavy-ion collisions (see e.g. [1]). It is characterized 
by the Fourier decomposition of the azimuthal particle 
distribution with respect to the reaction plane 



di<p 



\J ( 

= ^ 1 + 2 J] v„ cos [n(0 - >Prp)] 



, (1) 



where A^o is the number of particles, v„ is the n-th har- 
monic coefficient, (p is the azimuthal angle of particles 
and *Prp is the azimuthal angle of the true reaction plane 
which is determined by the beam axis and impact pa- 
rameter 

One of the standard methods to extract v„ is the event 
plane method 13]. The most important task in this 
method is to estimate the reaction plane from the mea- 
sured particles for each harmonic n. We do not dis- 
tinguish in this paper between reaction plane and par- 
ticipant plane since we discuss only resolution correc- 
tion. The estimated reaction plane is defined as the event 
plane ^„ (-n/n < ^'„ < n/n), but due to the finite mul- 
tiplicity in nuclear collisions, the event plane can be dif- 
ferent from the reaction plane. The observed v°''^(M) 
for a given, small, centrality range M must be corrected 
by the event plane resolution 'R„iM) in order to take into 
account the difiference between true reaction plane and 
event plane 
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«cos [n(0->PJ]))M 



(cos [n{^>„ - »Prp)]>m 

Preprint submitted to Physics Letter B 



KiM) ' 



(2) 



where m is the harmonic of the event plane and n = km 
is the harmonic of interest. Brackets denote the average 
over events, while double brackets denote the average 
over particles in all events. Subscript M on the bracket 
emphasizes that the average is taken for a given central- 
ity M. For simplicity, we will omit (M) from the observ- 
ables, for example, we will write A^o instead of Nq{M). 
In experiment the reaction plane angle Trp in Eq. (|2]i is 
not known. Therefore at least two subevent planes are 
necessary in order to calculate the event plane resolu- 
tion An average v„ over a wider centrality range 
R can be calculated once v"*"* and are determined 
within R 
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We introduced brackets (...) for simplicity, which denote 
the average over a wide centrality range weighted by 
particle multiplicity A'o for a given phase space (e.g. for 
a given transverse momentum). 

fin depends on multiplicity and v„ itself, therefore 
Eq. O) requires that vf' and are measured in suf- 
ficiently small centrality intervals. However, for rare 
particles (e.g. £1, J/t//) it is not always possible to do so. 
One of the conventional approaches (e.g. see Ref. f^) IS 
to average v"*"* as well as 'R„ separately in a wide central- 
ity range, weighted by the corresponding particle yields, 
and then take the ratio of ^v"*"*^ to (!R„). However, this 
approach systematically overestimates v„ as we will dis- 
cuss in the Section[3] The main point of this paper is the 
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following inequality 



where the left hand side is the correct average over a 
wide centrality bin, while the right hand side shows a 
commonly used approximation. 

We will show the proper way to correct for the finite 
event plane resolution in wide centrality bins. In Sec- 
tion |2l we derive the equations used to calculate (v'„) 
in our approach. We also show that the derived equa- 
tions are equivalent to the average calculated from nar- 
row centrality bins (see Eq. (|3]l). In Section|3] we show 
a fast Monte Carlo simulation to demonstrate the vahd- 
ity of the method, for the case of V2. 



2. Implementation 

We now show several practical implementations to 
correct v„ for the event plane resolution in wide cen- 
trality bins. There are two or three steps to calculate v„ 
in wide centrality bins 

1 . Determine the event plane resolution %, as a func- 
tion of narrow centrality ranges M. 

2. Analyse v„ with the weights 1 for any, wide, 
centrality range of interest. 

3. Perform an additional resolution correction (de- 
pends on implementation). 

In the following two subsections, we discuss detailed 
implementations how to apply corrections separately for 
different types of particle identification. For the sake of 
simplicity, we assume that non-flow effects are negligi- 
ble, and all correlations between particles are induced 
by flow. Namely, the azimuthal particle distribution 
with respect to the event plane can be written as 
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2.1. Event-by-event particle identification 

If the particle of interest can be identified on 
an event-by-event basis, one can directly calculate 
cos [n{(p - ^m)] for every particle, corrected with the 
event plane resolution for the corresponding centrality 
M 



cos [n{(f> - "¥,„)] 
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where is supposed to be averaged over many events 
in advance. The event and centrality average of term (|6) 
over <p - reduces in this case to Eq. ^ 
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The main difference of our implementation to the con- 
ventional approach is the event-by-event correction 
in Eq. O. As we have mentioned earlier, event-by- 
event correction does not mean that the should 
be calculated event-by-event, but rather the correction 
is made event-by-event. Practically, the statistical error 
on the event plane resolution can be ignored if the event 
plane resolution is good. However, if the event plane 
resolution is very small then one needs to propagate also 
the statistical error on the event plane resolution into the 
final Vn values. 

2.2. Statistical particle identification 

There are two approaches in case the particle yield of 
interest can only be extracted statistically: the invariant 
mass fit method and the event plane method. The in- 
variant mass fit method is almost equivalent to the one 
introduced in Section lZTl while the event plane method 
needs one additional step to obtain the final v„, which 
will be discussed below. 

2.2.1. Invariant mass fit method 

The invariant mass method is quite useful to analyse 
particles that are detected through their decay products, 
such as /Tj — > 71^ , A — > pn^ and so on. The point of 
this method is to calculate the mean cos [n(0 - '^„^} as 
a function of invariant mass Minv 
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where S is the signal yield, B is background yield, 
v^, and v^^* are the v„ for signal, background and 
total particles, respectively. Signal and background 
contributions are decomposed by taking into account 
the measured signal-to-background ratio and using a 
parametrization of the background shape |4]. Since 
the average cosine is calculated in this approach, one 
can directly extract the by subtracting the back- 
ground contribution. The only modification is to add 
a weight 1 /%, on an event-by-event basis when one 
fills the histograms for cos [n((f> - versus invariant 
mass similar to the Eq. Q. 

2.2.2. Event plane method 

The first step in the event plane method is the signal 
extraction for a given - bin 
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where N'^{<p - is the number of particles for a given 
<p - *P,„ bin, weighted for each centrality bin with the in- 
verse of the event plane resolution. The difference to the 
conventional method is the weight 1 /"Rn on the particle 
yields which will properly take into account the central- 
ity dependence of the event plane resolution. Second, 
one integrates Eq. (fTOl l over - *F,„ to calculate v„ 
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Note that v° ^' in Eq. ( fTTT i is different from the con- 
ventional v°''^ due to the !/■??„ weight on the particle 
yields. One could immediately notice that the numer- 
ator is identical to Eq. Q with the additional normal- 
ization 1 / (in dMNoj- The denominator in Eq. (fTTT i 
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Since the normalization factor 2n (j^ dMNo^ is can- 
celled out between numerator and denominator, the (v„) 



is obtained by multiplying (l/K) by the (v°'"''^^ 

(vr)-(vr'^)(^) = <v„). (13) 

We would Uke to emphasize two important points for 
the event plane method. First, one must take an average 
of the inverse event plane resolution 1 /%, (not an aver- 
age of Kn)- Second, the integral dM in Eq. ( fTTT i for 
both numerator and denominator should be calculated 
in the same phase space. For instance, if one measures 
the v„ as a function of transverse momentum pr, one 
should calculate {!/%,} for each pr bin. 

The new implementation can also be applied for the 
scalar product method IH 0] because the scalar prod- 
uct method is a simple extension of event plane method 
where the magnitude of the Q-vector is considered for 
both vf*^ andK- 



3. Simulation results 

In this section we validate our implementation by a 
fast Monte Carlo simulation for viipj)- As we have 
mentioned earlier, this approach can be applied for any 
harmonic of interest. 
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Figure 1 : Input multiplicity distribution in the fast Monte Carlo sim- 
ulation. The shaded area show 0-5%, 5-10%, 75-80% centrality 
classes. 

Figure [T| shows the input multiplicity distribution in 
0-80% centrality for the 50 million generated events. 
The event centrality classes are determined by using the 
multiplicity distribution as it is typically done in the ex- 
perimental data. We divide the multiplicity distribution 
into 16 bins with 5% increments within 0-80%. 
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For every event a number of tracks according to the 
input multiplicity distribution was sampled. A Boltz- 
mann like pj distribution from 0.25-3 GeV/c was used. 
The input v^Cpr) for every centrality bin was fixed to a 
parametrization of the following form 



V^ipr)^Ail-e-n[ ^_^^l^_,,, -d). (14) 

The parameter A was increased towards more periph- 
eral centralities. Furthermore we have added for every 
event a Gaussian smearing on v^^ with a width of 0.05. 
Figure 12] depicts the used v™{pt) for the 16 centrality 
bins. The event Q-vector was calculated according to 
Ref. H. 
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Figure 2: (Color online) Input V2{pt) distributions for different cen- 
tralities. Each line con'esponds to a 5% centrality range, v^ipr) in- 
creases towards more peripheral events. The applied Gaussian smear- 
ing is not shown. 



Figure[3]shows the random subevent plane resolution 
as a function of centrality. For each event, particles are 
randomly divided into two different groups in order to 
evaluate the resolution The result shown here is the 
subevent plane resolution by using Eq. (14) in Ref. 
The resolution reaches a maximum of 0.45 around 30% 
centrality, it decreases towards more central and periph- 
eral events because of a lower V2 in central and less mul- 
tiplicity in peripheral events. 

Figure |4] shows an example of a particle azimuthal 
distribution with respect to the second harmonic event 
plane in a nan^ow pj bin. The particle yields are 
weighted by the inverse event plane resolution as indi- 
cated by in the y-axis title. The dashed line rep- 
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Figure 3: (Color online) Random subevent plane resolution for the 
second harmonic plane as a function of centrality. The resolution is 
calculated from random subevents. 
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Figure 4: (Color online) Azimuthal distribution of particles with re- 
spect to the second harmonic event plane at 1.5 < pj < 1.56 GeV/c 
in 0-80% centrality bin. The dashed line is the fit result for 
l+2vf"cos [20-2'F2]. 
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resents a fit with 1 + Ivf""'"^ cos [20 - I'Vil The yield 
extraction and fit are repeated for all other pj bins. 

Figure |5] shows (V2) as a function of pr in the 0-80% 
centrality bin. For comparison, we also plot the ob- 
served ^v*'*^ without resolution correction as shown by 
solid grey circles. The difference between observed and 
corrected (V2) gives the size of the resolution correction. 
The input ^v'"^ of the simulation is shown as a magenta 
solid line. We tested both, the direct cosine calculation 
^v™**^ from Eq. (jTjl and the event plane method (v'^, '^'^ 
from Eq. (fT3T l. as shown by open blue stars and solid 
blue stars, respectively in panel (a) and (b). One can 
see that both direct cosine and event plane methods are 
consistent with the input (v"). The data points are in 
agreement with within 0.3%. Panel (b) shows the ob- 
served (vj**^^, corrected by the inverse mean event plane 
resolution ' (green line) and the mean inverse event 
plane resolution ('^ ') (red line). In panel (c) we show 
a correction similar to panel (b) but with the event plane 
resolution calculated for every single pj bin. The rela- 
tive difference A (V2) of the corrected (V2) values to the 
input ^v™^ is shown in the lower panels. For the correc- 
tion with (^)"' and we get a A<V2) / (v^") of 5% 
and 15%, respectively. The difference becomes smaller 
if the average resolution is calculated for every pj bin 
as shown in panel (c), but we still obtain relative devi- 
ations of 4.5% {{K)-]:) and 10% ((•;?"') ). In all cases 
the relative difference is almost independent of pj. 

4. Summary 

We have introduced an implementation to avoid event 
plane resolution effects of wide centrality bins on the 
v„ measurements. The new approach is essentially an 
extension of existing event plane methods. The basic 
modification is to add an inverse of event plane resolu- 
tion weight to the cosine term or to the particle yields 
depending on the method. We confirmed that our ap- 
proach reproduces the input ^V2'^(/?7-) in a wide central- 
ity bin, whereas the conventional corrections lead to a 
bias as large as -15% on the reconstructed (v2> with 
our fast Monte Carlo simulation. 
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Figure 5: (Color online) Results of V2 as a function of pr with differ- 
ent approaches of resolution corrections. The magenta solid line cor- 
responds to the input (V2), gray filled circles are the observed ^i'"*"^^ 
without resolution connection, blue open stars in panel (a) are from 
the direct cosine calculation from Eq. (??), blue filled stars in panel 
(b) are from the event plane method from Eq. ( 1131 . Other lines show 
the observed (vf'^^ with different mean event plane resolution correc- 
tions. In panel (c) the event plane resolution is calculated for every px 
bin separately. 
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